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I.  INTRODUCTION 


The  flight  of  projectiles  covers  a  wide  range  of  speeds.  The  capability 
to  predict  the  aerodynamic  behavior  in  transonic  and  supersonic  flow  is  essen¬ 
tial  to  the  design  of  new  projectiles  or  in  the  modification  of  existing 
shapes.  In  recent  years  a  considerable  research  effort  has  been  focused  on 
the  development  of  modern  predictive  capabilities  for  determining  projectile 
aerodynamics.  The  general  approach  taken  has  been  to  apply  the  thin-layer 
Navier-Stokes  solvers  in  both  a  space  marching  and  time  marching  manner  to  the 
supersonic  and  transonic  flow  fields  respectively. 

At  supersonic  speeds  Sturek  and  Schiff^  have  successfully  applied  a 
Parabolized  Navier-Stokes  code  to  predict  the  flow  field  about  an  ogive- 
cylinder-boattail  body.  The  resultant  aerodynamic  coefficients  were  found  to 
be  in  very  good  agreement  with  the  available  experimental  data.  This  work  has 
been  further  advanced  to  include  the  effects  of  a  blunt  nose  configuration^ 
which  more  closely  resembles  an  actual  fielded  artillery  round.  The  initial 
emphasis  has  been  on  developing  a  specific  capability  for  predicting  projec¬ 
tile  aerodynamics.  This  computational  capability  is  most  useful  in  the  even¬ 
tual  design  and  evaluation  process.  Towards  this  end  a  computational  param¬ 
etric  study  has  been  performed  by  Sturek  and  Mylin^  and  trends  in  the  aerody¬ 
namic  coefficients  were  determined. 

The  transonic  regime  represents  an  area  in  which  a  similar  computational 
capability  is  being  sought.  The  critical  aerodynamic  behavior  which  occurs 
between  Mach  =  .90  and  1.2  is  always  of  concern  in  the  development  of  shell. 
For  the  transonic  regime  the  time  dependent  thin-layer  Navier-Stokes  technique 
has  been  applied  to  projectile  configurations  for  axisymmetric  flow  condi¬ 
tions. A  base  flow  version  of  the  code  has  also  been  developed. 5  Agreement 
with  experimental  data  for  the  surface  pressures  has  been  very  good.  The  base 
flow  calculations  shcx'v  correct  qualitative  features  of  the  flow  field.  The 
extension  of  these  techniques  to  a  fully  three  dimensional  flow  field  which 
exists  on  a  spinning  projectile  at  angle  of  attack  has  recently  been 
reported.^  The  normal  force  coefficient  was  Shown  to  be  in  good  quantitative 
agreement  with  the  available  experimental  data.  The  Magnus  force  predictions 
were  not  as  encouraging  but  did  show  the  correct  development  of  Magnus  force 
along  the  body. 

This  computational  study  utilizes  the  developed  time  dependent  thin-layer 
Navier-Stokes  computational  technique  to  compute  the  static  and  dynamic  aero¬ 
dynamic  coefficients  for  a  standard  M549  artillery  shell.  The  initial  objec¬ 
tive  was  to  obtain  a  series  of  computational  results  which  cover  the  critical 
Mach  number  range  and  determine  the  applicability  of  these  codes  as  a  predic¬ 
tive  tool.  Three  versions  of  the  thin-layer  Navier-Stokes  codes  were  used: 
(a)  a  full  three  dimensional  version  (F3D);  (b)  vectorized  symmetry  version 
(VSYMSD),  and  (c)  a  axisymmetric  base  flow  version.  The  full  range  of  aerody¬ 
namic  coefficients,  which  includes  base  drag  have  been  determined.  The  com¬ 
puted  results  are  compared  to  a  set  of  experimental  data  which  combines  both 
wind  tunnel  and  range  data.  All  computations  for  which  the  aerodynamic  coef¬ 
ficients  are  presented  were  obtained  on  the  Sandia  National  Laboratory  Cray  IS 
which  has  one  million  words  of  memory.  Results  from  computations  obtained  on 
the  2H  word  Cray  IS  at  Lawrence  Livermore  National  Laboratory,  are  also  shown. 
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II.  GOVERNING  EQUATIONS  AND  COMPUTATIONAL  TECHNIQUE 

The  governing  equations  solved  are  the  unsteady  thin-layer  Navler-Stokes 
equations.  The  thin-layer  approximation,  which  neglects  the  viscous  gradients 
in  the  longitudinal  and  circumferential  direction,  has  been  applied  to  a 
variety  of  problems  with  very  good  results.  The  transformed  equations  written 
in  strong  conservation  law  form  for  the  full  3D  (Eq.  (1))  and  axi symmetric 
formulation  (Eq.  (2))  are 


aq  .  3E  .  3F  .  3G  _  1  3S 
TF  Tl  *50  Tc  “  Re  3c 
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respectively.  The  general  coordinate  transformations  are  defined  as 


s 
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^(x.y.z.t) 

n(x,y,z,t) 

c(x,y,z,t) 


=  is  the  longitudinal  coordinate 
=  is  the  circumferential  coordinate 
=  is  the  near  normal  coordinate 


T  =  time 


The  vector  q  contains  the  dependent  variables  (p,  pu,  pv,  pw,  e)  and  the  flux 

AAA 

vectors  E,  F,  G  contain  terms  which  arise  from  the  conservation  of  mass, 

A 

momentum  and  energy  in  the  three  coordinate  directions.  The  source  vector  H 
contain  terms  which  result  from  an  analytic  determination  of  the  flux  vector 

F  given  the  assumptions  of  axisymmetric  flow  and  a  constant  angular  velocity. 

A 

The  viscous  terms  are  contained  in  the  vector  S  which  is  seen  to  have  varia¬ 
tion  in  the  c  direction  only.  This  represents  the  thin-layer  approximation 
mentioned  earlier. 

The  numerical  scheme  used  for  the  solution  of  Equations  (1)  and  (2)  is  a 
fully  implicit,  approximately  factored,  finite  difference  algorithm  in  delta 
form  as  analysed  by  Beam  and  Warming.^  The  solution  is  implemented  by  an 
approximate  factorization  which  allows  the  system  of  equations  to  be  solved  in 
three  successive  one  dimensional  steps  for  the  3D  problem  and  two  one  dimen¬ 
sional  steps  for  the  axisymmetric  formulation.  DetaiOs  of  the  numerical 
method,  computational  algorithm  and  boundary  conditions  can  be  found  in 
References  6  and  8. 
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III.  MOn^L  AND  COMPUTATIONAL  GRIDS 

t 

The  actual  projectile  configuration  for  which  these  simulations  were  run 
IS  shown  in  Figure  1.  There  are  certain  features  on  this  projectile,  such  as 
the  flat  nose  and  rotating  band,  which  can  not  yet  be  modeled  exactly.  There¬ 
fore  a  modified  configuration  as  shown  In  Figure  2,  was  used  for  the  computa¬ 
tions.  The  flat  nose  was  modeled  as  a  hemisphere  cap  and  the  rotating  band 
was  omitted.  Recent  results  have  been  presented  by  Danberg^  for  supersonic 
flow  wherein  the  rotating  band  was  modeled  as  a  ramp  function  and  the  results 
compared  well  the  experimental  data.  This  work  Is  presently  being  extended  to 
for  the  transonic  regime. 

Equations  (1)  and  (2)  are  solved  In  a  transformed  coordinate  space  thus 
allowing  for  a  wide  variation  in  body  geometry.  The  body  configuration  is 
introduced  through  the  determination  of  the  metric  terms  etc.  and 

the  Jacobian.  These  terms  are  formed  by  a  combination  of  the  derivative  terms 
x^,  y^,  2^.  One  of  the  first  steps  in  performing  a  computation  is  the  genera¬ 
tion  of  a  stationary  computational  grid  for  the  metric  determination.  These 
points  are  determined  prior  to  the  computations  and  are  not  changed  with  time. 
Three  grids  were  used  in  this  study  for  the  3D  calculations  and  one  grid  for 
the  base  flow  results.  All  30  grids  were  generated  using  a  projectile  grid 
generation  code^°  which  solves  either  elliptic  or  by  hyperbolic  differential 
equations. 

The  grid  used  for  the  full  3D  run  is  shown  in  Figures  3  and  4  consisted 
of  60  longitudinal,  24  normal  and  36  circumferential  points.  The  overall  grid 
was  extended  to  18  calibers  in  front,  above  and  behind  the  projectile  body. 
This  is  required  to  insure  that  imposed  free  stream  boundary  conditions  are 
satisfied  and  that  waves  are  not  reflected  back  into  the  computational  domain. 
This  type  of  computational  space  has  also  been  verified  by  a  series  of  axlsym- 
metric  calculations.  An  expanded  view  of  the  grid  is  shown  In  Figure  4  and 
shows  the  details  of  longitudinal  and  radial  grid  point  clustering  required 
for  the  shock  and  viscous  resolution  respectively.  The  wake  has  been  modeled 
as  an  extended  sting  for  the  3D  calculations.  The  base  flow  code  models  the 
base  regime  and  calculates  the  wake  and  resulting  base  pressure, 

A  second  grid,  shown  in  Figure  5,  was  developed  for  the  symmetry  calcula¬ 
tions  and  contained  80  longitudinal,  32  normal  and  21  grid  points  in  the 
circumferential  half  plane.  Only  the  static  aerodynamic  coefficients  can  be 
determined  with  this  grid  since  bi-lateral  symmetry  conditions  are  Imposed. 
Thus  the  full  three  dimensional  flow  field  which  exists  for  a  spinning  projec¬ 
tile  at  angle  of  attack  cannot  be  determined.  Grid  points  for  both  cases  were 
clustered  on  the  boattail  to  resolve  the  expected  asymmetric  shock  pattern. 
The  shock  asymmetry  and  possible  base  region  influence  are  felt  to  be  the 
major  contributions  to  the  critical  aerodynamic  behavior.  Only  the  static 
aerodynamic  coefficients  can  be  determined,  however  since  this  grid  Imposes  a 
bi-lateral  symmetry  condition.  Thus  the  full  three  dimensional  flow  field 
which  exists  for  a  spinning  projectile  at  angle  of  attack  cannot  be  deter¬ 
mined.  These  two  grids  made  use  of  the  full,  one  million  word  memory  capacity 
on  the  Sandia  Cray  IS. 

A  third  grid  was  generated,  for  a  symmetry  calculation  run  on  the 
Lawrence  Livermore  Cray  IS  with  two  million  words,  and  Is  shown  in  Figure  6. 
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This  machine  enabled  a  grid  of  110  longitudinal,  32  normal  and  21  circumfer¬ 
ential  points  to  be  generated.  The  fine  resolution  in  the  longitudinal  direc¬ 
tion  is  evident  from  this  figure  and  improved  results  are  expected.  The  aero¬ 
dynamic  coefficients  for  this  last  case  have  not  yet  been  computed  however,  a 

preliminary  result  of  the  pressure  distribution  along  the  surface  will  be 
shown. 

A  final  grid,  generated  using  a  hybrid  grid  generation  code,  was  used  for 
the  base  flow  calculations  and  is  shown  in  Figure  7.  Grid  points  are  shown  to 
be  clustered  in  the  base  region  and  although  it  has  not  been  used  here  the 
capability  exists  to  have  the  base  grid  adapt  to  the  shear  layer  as  it 
develops. 


IV.  COMPUTATIONAL  CONDITIONS  AND  CODES 

An  objective  of  this  effort  was  to  computationally  determine  the  magni¬ 
tude  and  trends  of  the  static  and  dynamic  aerodynamic  coefficients  where  the 
critical  behavior  is  known  to  exist,  namely  between  Mach  =  .90  through  1.2. 
The  following  table  is  a  summary  of  the  conditions  and  codes  used  for  the 
computati ons. 


MACH  NO. 

(Re  X  10"6) 

F3D 

VSYM3D 

BASE  FLOW 

.90 

2.2 

X 

X 

X 

.95 

2.32 

X 

X 

X 

• 

LO 

oo 

2.39 

X 

X 

X 

1.05 

2.56 

X 

X 

1.10 

2.69 

X 

X 

1.20 

2.93 

X 

V 

/V 

A  large  number  of  runs  were  made  during  this  study  using  the  full  3D  and 
vectorized  symmetry  30  Navier-Stokes  code  on  the  Sandia  National  Laboratory 
Cray  IS.  The  computations  for  base  pressure  were  all  performed  at  the 
Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland  on  the  CDC 
7600.  The  run  times  for  the  30  cases  were  quite  large.  The  F3D  code  ran  at 
.00055  sec/time  step/point  and  required  1000  steps  for  convergence.  Thus  each 
case  ran  about  6.9  hours  on  the  Cray.  The  VSYM3D  ran  at  .00011  sec/time  step/ 
point  and  required  2500  time  step  for  convergence.  The  total  run  time  per 
case  for  this  version  was  3.1  hours. 

The  two  codes  differed  in  the  operational  aspect  with  the  F3D  being  more 
robust  than  the  VSYM30.  Lower  values  of  smoothing  parameters  (explicit 
smoothing/At  ■  1.6)  and  higher  values  of  At  {•*  .05)  were^able  to  be  used  for 
the  F3D  code  while  a  smoothing  value  of  -  4.0  and  maximum  At  •  .03  were 
required  for  the  VSYM3D  code.  The  lower  At  thus  caused  the  VSYM3D  to  require 
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more  time  steps  for  convergence.  Additlonally  the" cohvefgenfe  Histories’ for 
the  F3D  runs  were  favorable  in  that  the  residual  decrease  and  time  Increases 
did  not  cause  the  solution  to  diverge  repeatedly.  A  very  careful  approach  to 
changing  times  and  smoothing  was  required  for  the  VSYM3D  runs  to  avoid  the 
solution  becoming  unstable.  However,  for  the  cases  with  Hach  >  1  the  F3D  code 
became  very  unstable  and  solutions  were  not  obtained.  The  cause  of  the  Insta¬ 
bility  was  due  to  the  shock  system  forming  in  front  of  the  nose.  Additional 
local  smoothing  would  have  been  required  to  obtain  a  solution. 


V.  RESULTS 

The  VSYM3D  code  was  used  to  determine  the  static  coefficients  only,  since 
bi-lateral  symmetry  was  imposed  and  the  projectile  could  not  spin.  The  F3D 
code  was  used  to  determine  the  flow  fields  for  the  spinning  cases  from  which 
the  static  and  dynamic  coefficients  could  be  obtained.  Qualitative  features, 
resulting  from  the  F3n  calculations,  are  shown  in  Figures  8a  and  8b,  where  the 
Mach  contours  have  been  plotted  for  M  =  .90  and  M  =  .95  respectively.  The 
supersonic  regions  are  clearly  seen  in  these  figures  as  well  a.s  the  asymmetry 
of  the  shock  positions  due  to  the  projectile  being  at  angle  of  attack.  The 
longitudinal  asymmetry  and  shock  interaction  with  the  base  region  flow  field 
are  believed  to  cause  the  critical  behavior  in  the  aerodynamic  coefficients 
over  the  transonic  range. 

The  resultant  flow  field  calculations  were  integrated  to  determine  the 
aerodynamic  coefficients.  The  results  to  be  presented  here  include  the  static 
coefficients:  normal  force  slope,  pitching  moment  slope,  and  center  of  pres¬ 
sure.  The  dynamic  coefficients  presented  include:  Magnus  moment  slope,  side 
force  center  of  pressure,  and  roll  moment  slope.  The  sign  convention  used  is 
shov/n  in  Figure  9.  The  normal  force  is  shown  to  act  in  the  plane  of  the  angle 
of  attack  while  the  Magnus  force  acts  in  the  plane  normal  to  the  angle  of 
attack.  The  Magnus  moment  is  shown  to  act  about  the  normal  force  axis  while 

the  pitching  moment  acts  about  the  side  force  axis.  The  slope  of  the  moment 

coefficients  (C^,.  ,  )  are  presented,  however,  since  it  is  those  terms  which 

pa  a 

are  used  to  dete'-mine  stability  criteria.  The  computed  results  are  compared 
with  experimental  data^2  which  are  a  compilation  of  wind  tunnel  and  range  data 
for  the  XMd49  shape. 

The  static  aerodynamic  coefficients  for  normal  force  slope,  pitching 
moment  slope  and  center  of  pressure  are  shown  as  a  function  of  Mach  number  in 
Figures  10,  11,  and  12  respectively.  The  results  from  both  codes,  the  F3D  and 
VSYM3D  are  compared  with  the  experimental  data.  The  primary  difference  in  the 
two  codes  to  be  considered  here  is  the  total  number  and  distribution  of  grid 

points.  The  VSYM3D  version  allows  more  points  to  be  used  in  the  longitudinal 

and  normal  directions  than  that  for  the  F3D  version.  The  circumferential 
point  distribution  has  been  fixed  at  =  10®  for  both  versions.  The  normal 
force  slope  has  been  shown  to  depend  on  the  number  of  circumferential  points 
by  Klopfer  and  Chaussee. T^-js  work  showed  that  the  normal  force  slope 
continued  to  change  until  points  were  used  every  10“  in  the  circumferential 
plane.  That  same  criteria  was  used  for  this  work,  ^ 

The  normal  force  slope  coefficient  is  shown  in  Figure  10  for  both  the 
computational  results  and  experimental  data.  The  experimental  data  is  shown 
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to  decrease  as  it  approaches  a.  critical  Mach  .'yrher  tM  -and  then  shows  a 
rapid  increase  throughout  the '.remaining  transonic  regime.  The  F3D  computa¬ 
tional  results  agree  well  above  M  =  .95  however,  it  over  predicts  the  data 
below  M  =  .95.  The  VSYM3D  results  show  an  improved  agreement  with  the  experi¬ 
mental  Gita  except  at  M  =  1.2.  The  increased  number  of  grid  points  was 
expected  to  provide  better  agreement  and  generally  has  predicted  the  correct 
trend  for  C|^  . 

a 

The  pitching  moment  slope  coefficient  is  gerterally  of  greater  concern  in 
projectile  aerodynamics.  Figure  11  shows  the  pitching  moment  slope  comparison 
for  the  computations  and  the  experimental  data.  Here  the  F3D  results  are 
found  to  be  in  general  agreement  with  the  magnitude  of  the  experimental  data 
but  show  very  little  effect  of  Mach  number.  The  inability  to  predict  the  Mach 
number  effect  is  attributed  to  the  insufficient  number  of  grid  points  along 
the  body.  The  VSYM3D  computations  howeve’*  show  a  very  definite  trend  with 
Mach  number.  This  is  a  result  of  being  able  to  capture  more  accurately  the 
expansion  regions  and  shock  structure  with  the  increased  resolution.  Although 
the  magnitude  of  the  VSYM3D  results  do  not  agree  with  the  experimental  data, 
the  prediction  of  the  trend  with  Mach  number  is  very  encouraging  and  was  one 
of  the  objectives  of  this  study. 

The  total  drag  of  a  projectile  configuration  results  from  summation  of 
the  pressure,  viscous,  and  base  drag  components.  The  30  codes  were  used  to 
obtain  the  pressure  and  viscous  contributions  to  drag  since  the  base  was 
modeled  as  an  extended  sting.  The  base  drag  was  obtained  by  running  the  base 
flow  (axi symmetri c)  code  at  a  =  0.0“.  The  total  drag  was  then  determined  by 
adding  all  the  components  described  above.  The  only  effect  not  accounted  for 
in  this  approach  is  the  change  in  base  drag  due  to  angle  of  attack.  Figure  12 
compares  the  total  drag  coefficients  for  both  computations  against  the  experi¬ 
mental  data.  The  F30  results  do  show  a  reasonable  agreement  with  the  data 
however  they  are  inconclusive  due  to  the  limited  number  of  calculations  per¬ 
formed  in  this  critical  Mach  number  region.  The  VSYM3D  results  cover  a  wider 
range  and  tends  to  predict  the  drag  rise  associated  with  transonic  flow.  The 
magnitude  at  K  =  1.05  and  M  =  1.2  does  not  agree  with  the  experimental  data. 
Although  the  pressure  drag  component  had  increased  at  M  =  1.2,  the  base  drag 
was  still  predicting  a  significant  contribution.  The  difference  between  the 
experimental  data  and  calculations  are  most  likely  due  to  the  overprediction 
of  the  base  drag. 

The  dynamic  coefficients  of  Magnus  moment  slope,  side  force  center  of 
pressure  and  roll  moment  are  shown  in  Figures  13,  14,  and  15  respectively. 
The  flow  field  computations  for  these  coefficients  required  the  F3D  code  since 
a  projectile  at  angle  of  attack  with  spin  is  a  fully  three  dimensional  problem 
with  no  plane  of  symmetry.  The  critical  aerodynamic  behavior  is  shown  experi¬ 
mentally  to  exist  for  the  dynamic  coefficients  as  well.  The  experimental 
results  presented  in  Figure  13  show  the  Magnus  moment  slope  to  change  signifi¬ 
cantly  throughout  the  transonic  regime.  The  computed  results  however  show 
almost  no  change  in  the  Magnus  moment  coefficient  with  Mach  number.  This 
large  change  in  aerodynamic  coefficients  is  in  part  attributed  to  the  longitu¬ 
dinal  asymmetry  of  the  shock  position.  This  has  not  been  captured  by  the 
present  computation  due  to  the  poor  grid  resolution  in.,the  boattail  region. 
The  Magnus  force  center  of  pressure  is  shown  in  Figure  14  and  although  shows 
some  change  with  respect  to  the  Mach  number  variation  is  considered  to  be  only 
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a  marginal  result.  The  Magnus  effect  is  a  result  of  the  viscous  interation 
and  resulting  skewed  displacement  surface  which  exists  about  a  spinning  pro¬ 
jectile  at  angle  of  attack.  Therefore  adequate  resolution  of  the  viscous 
boundary  layer  is  essential  for  good  results.  An  indication  of  the  inadequate 
viscous  resolution  for  these  computations  is  shown  in  Figure  15  where  the  roll 
moment  coefficient  has  been  plotted  as  function  of  Mach  number.  The  roll 
moment  is  strictly  a  result  of  the  viscous  effects  and  thus  is  a  much  more 
sensitive  indicator  of  the  viscous  resolution  than  the  Magnus  force.  The  poor 
agreement  shown  here  is  therefore  directly  attributed  to  the  insufficient  grid 
resolution  in  the  boundary  layer.  Using  the  full  capacity  of  the  CRAY  IS 
allowed  for  only  24  points  between  the  body  and  the  outer  boundary  of  18 
calibers. 

The  longitudinal  distribution  of  points  is  also  very  critical  to  the 
determination  of  shock  position  and  eventually  to  the  determination  of  the 
pressure  acting  on  the  body  surface.  Figure  16  is  a  plot  of  the  accumulative 
pitching  moment  as  a  function  of  position  on  the  model  for  various  Mach 
numbers.  The  pitching  moment  on  the  nose  is  seen  to  be  the  same  for  all  Mach 
numbers.  The  largest  effect  is  shown  to  occur  at  the  ogive  cylinder  junction 
and  on  the  boattail.  With  the  boattail  having  such  a  dramatic  effect  on  the 
aerodynamic  coefficients,  good  resolution  of  the  boattail  flow  field  is  essen¬ 
tial.  A  few  runs  were  made  on  the  Lawrence  Livermore  Laboratory,  CRAY  IS  with 
two  million  words,  using  the  grid  shown  in  Figure  6.  A  plot  of  the  pressure 
distribution  acting  on  the  leeward  surface  is  presented  in  Figure  17  for  the 
F30,  VSYH3D  (one  million)  and  VSM30  (two  million)  runs.  With  the  increased 
resolution  a  sharper  definition  of  the  shock,  indicated  by  a  more  rapid  rise 
in  pressure,  is  seen  on  the  cylinder  (X/0  «  3.8)  and  on  the  boattail 
(X/D  »  5.5).  This  pressure  distribution  is  expected  to  provide  a  more 
accurate  prediction  of  the  aerodynamic  coefficients. 


VI.  CONCLUDING  REMARKS 

Three  versions  of  the  unsteady  thin-layer  Navier-Stokes  codes  were 
applied  to  the  standard  M549  projectile  configuration.  The  full  3D  version, 
F3D,  was  used  to  determine  the  static  and  dynamic  aerodynamic  coefficients  as 
a  function  of  Mach  number.  The  vectorized  symmetry  3D  version,  VSYM3D,  was 
used  to  determine  the  static  coefficients.  A  base  flow  version  was  used  to 
compute  the  base  drag.  The  computed  coefficients  have  been  compared  to  a  set 
of  experimental  data.  The  static  coefficients  predicted  the  correct  trends 
with  the  experimental  data  while  the  dynamic  coefficients  did  not.  The 
critical  aerodynamic  behavior  was  found  to  be  qualitatively  predicted  when 
sufficient  longitudinal  grid  points  were  made  available  by  using  the  symmetry 
version  of  the  code.  Additional  work  is  required  in  resolving  the  viscous 
effects  of  a  spinning  projectile  at  angle  of  attack. 
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Figure  12.  Total  drag  coefficient  vs  Mach  number. 
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Figure  13.  Wagnus  moment  slope  coefficient  vs  Mach  number 
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Figure  lA.  Side  force  center  of  pressure  vs  Mach  number. 
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Figure  15.  Roll  moment  coefficient  vs  Mach  number 
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Figure  P.  Leeward  surface  pressure  vs  longitudinal  position,  F30 
_ VSYM3D(iM).  VSYM3D(2M),  M  =  .95.  a  =  2.0. _ 
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LIST  OF  SYMBOLS 


a 


CP 

CPY 


d 


0 

A  A  ^ 

E,F,G 

H 


L 

M 


P 

A 

q 

anTv 


Re 


S 


u,v,w 


=  speed  of  sound 

=  center  of  pressure  for  normal  force 

=  center  of  pressure  for  Magnus  force 

=  roll  damping 

=  pitching  moment  coefficient 

=  d2c^/[d(-^)»da],  slope  of  Magnus  moment  coefficient  evaluated  at 

PD/V  =  n,  a  =  0 

=  dC^/da,  slope  of  pitching  moment  coefficient  evaluated  at  a  =  0 

=  Magnus  (yawing)  moment  coefficient 

=  normal  force  coefficient 

=  dC,^/dci,  slope  of  normal  force  coefficient  evaluated  at  a  =  0 

=  Magnus  (side)  force 
=  local  diameter  of  model 
=  reference  diameter  of  model 

=  flux  vectors  of  transformed  gasdynamic  equation 
=  source  term  vector 
=  reference  length 
=  Mach  number 
=  pressure  normalized  by 
=  vector  of  dependent  variables 
=  non-dimensional  spin  rate  about  model  axis 
=  Reynolds  number 
=  viscous  flux  vector 

=  Cartesian  velocity  components  along  the  x,  y,  z  axis,  respectively, 
normalized  by  a„ 
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LIST  OF  SYMBOLS  (Continued) 


angle  of  attack 


C.n,; 


computational  coordinates  in  the  axial,  circumferential,  and  radial 
di recti ons 


=  density,  normalized  by  free-stream  density  p. 


Subscri pts 


free-stream  conditions 

spin  rate  about  model  axis,  rad/sec 
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